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Abstract: Some years ago, a method was proposed for measuring the CP-violating 
phase 7 using pairs of two-body decays that are related by U-spin reflection {d -H- s). 
In this paper we adapt this method to charmless B —)■ PPP decays. Time-dependent 
Dalitz-plot analyses of these three-body decays are required for the measurement of 
the mixing-induced CP asymmetries. However, isobar analyses of the decay ampli¬ 
tudes are not necessary. A potential advantage of using three-body decays is that the 
effects of U-spin breaking may be reduced by averaging over the Dalitz plot. This 
can be tested independently using the measurements of direct CP asymmetries and 
branching ratios in three-body charged B decays. 
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1. Introduction 

The standard model (SM) explanation of CP violation is that it is due to a phase in 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. This phase informa¬ 
tion is elegantly encoded in the unitarity triangle, whose interior CP-violating angles 
are a, [3 and 7 [1]. Using B decays, a great deal of effort has gone into measuring 
these angles in many different ways, along with the sides of the unitarity triangle, 
to search for inconsistencies that would indicate the presence of new physics (NP). 
Unfortunately, to date no such indications have been seen. This suggests that the NP 
is more massive than hoped for (which is consistent with the absence of NP signals 
at the LHC), and that the observation of its effects on CP violation in the B system 
will require measurements of greater precision. 

One interesting procedure for searching for NP involves the CKM phase 7 . The 
conventional way of measuring 7 uses the tree-level decay B~ —)■ [2-6]. Its 

latest value is 7 = (71.7ll;l)° [7]. However, suppose that 7 could be measured using 
decays that have significant (gluonic or electroweak) penguin contributions. If NP 
is present, it is likely to affect the (loop-level) penguins, in which case the extracted 
value of 7 would be different from that found using B~ —)■ . That is, one 

can probe NP by comparing the “tree-level” and “loop-level” values of 7 . (But note 
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that, if there is NP, the “loop-level” value of 7 will not be constant. It will generally 
vary, depending on which decays are used for its extraction.) 

One example of this involves B —)■ t^K decays. (In what follows, we briefly 
describe the method, but we refer the reader to Ref. [ 8 ] for full details.) There 
are four such decays: 7 r+iP°, B^ —)■ —)■ and B^ —)■ 

Using these processes, one can measnre nine observables: four branching ratios, fonr 
direct CP asymmetries, and one indirect (mixing-induced) CP asymmetry. However, 
assuming flavor SU(3) symmetry, the amplitudes can be written in terms of eight 
theoretical parameters: the magnitndes of the diagrams P^, T', C", three relative 
strong phases, and the weak phase 7 . (The valne of the weak phase (3 is taken from 
the measurement of indirect CP violation in B^{t) J/'^Ks [7].) With more 
observables than theoretical parameters, one can perform a £t to extract 7 . The 
valne fonnd is 7 = (35.3 ± 7.1)° [ 8 ], which differs from the tree-level value of 7 by 
3 . 5 ( 7 . While this is intriguing, one must remember that there is also an nnknown 
theoretical uncertainty due to SU(3) breaking. Before any conclusions can be drawn, 
there must be other, independent determinations of loop-level values of 7 . 

In 1999, R. Fleischer proposed a method for extracting 7 from P° —)■ K^K~ 
and —)■ 7 r’'' 7 r“, two decays whose amplitndes are related by U-spin {d -H- s) sym¬ 
metry [9]. Since penguin contributions are important for snch decays, this method 
would determine a loop-level valne of 7 . It reqnires the measnrement of the branching 
ratios and CP asymmetries, both direct and indirect, of both decays. This method 
is unaffected by final-state interactions; its theoretical accuracy is limited only by 
the size of U-spin-breaking effects. The factorizable U-spin-breaking corrections are 
calculable theoretically in terms of form factors and decay constants [9-11]. How¬ 
ever, the precise valne of the nonfactorizable U-spin-breaking correction is nnknown, 
thongh it may be sizeable [ 12 ]. 

Recently, the direct and indirect CP asymmetries in B^ —)■ K^K~ were measured 
by the LHCb Collaboration [13], and they carried out the above extraction of 7 [14]. 
Allowing for a U-spin-breaking error of 50%, they find 7 = (63.5lg7)°. However, if 
the theoretical error is > 60%, the nncertainty on 7 is mnch larger. 

It was pointed out in Ref. [9] that, with an additional dynamical assumption, one 
could replace P° —)■ K^K~ with B^ —)■ and analyses with this second decay 

were carried out in Refs. [10,11,15]. However, Ref. [16] hnds that the experimental 
data suggest that there may be a large nonfactorizable U-spin-breaking correction 
between B^ —)■ and B^ —)■ 7 r’'' 7 r“. This would lead to a large (unknown) 

theoretical error in the extraction of 7 using —)■ and —)■ 7 r''' 7 r“. 

The main purpose of the present paper is to note that the method of Ref. [9] can 
also be applied to charmless B PPP decays (P is a pseudoscalar meson) whose 
amplitudes are related by U spin. The key point is that, by using the Dalitz plots of 
the three-body decays, the effect of U-spin breaking may be greatly reduced. If this 
is possible - and there is an independent test to see if the procedure works - then 
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the loop-level value of 7 can be determined with little theoretical error. This will 
then provide a clean test for NP. 

Note that, nnder flavor SU(3) symmetry, the three hnal-state particles in charm¬ 
less B —)■ PPP decays are treated as identical, so that the six permntations of these 
particles must be considered. There have been a nnmber of papers recently that 
use the fully-symmetric hnal state [17-21], which can be obtained using an isobar 
analysis of the Dalitz plot. However, we stress that snch an analysis is not needed 
for the above method of extracting 7 - the fnll Dalitz plot is used. 

Examples of pairs of decays to which this method can be applied are (i) —)■ 

Kstt+tt- {b d) and 5° ^ KsK+R- {b ^ s), and (ii) 5° ^ KsK+R- {b d) 

and 5° ^ (b ^ s). The time-dependent Dalitz plots for B^ —)■ RsR^R~ 

and B^ Rs7r~^'n'~ were measured by the BaBar and Belle Collaborations [22-25], 
and a study of B^^-^ Rsh'^h'~ was made by the LHCb Collaboration [26]. For the 
B^g decays, it appears that B^g —)■ Rst^^'k~ is more promising experimentally. The 
hrst observation of this decay was reported in Ref. [26], and a stndy of the fntnre 
prospects for the measnrement of its time-dependent Dalitz plot was presented in 
Ref. [27]. Hopefnlly the method will be applied to decays B^ —)■ and 

—)■ RsR^R~ to extract 7 . 

In Sec. 2, we briefly discuss Dalitz plots and the distinction between the hnal 
states / and /. The U-spin relation between b ^ d and b ^ s decays is discussed in 
Sec. 3. In Sec. 4, we present the method for extracting 7 from a Dalitz-plot analysis 
of three-body decays. The subject of U-spin-breaking effects - the theoretical idea 
of how they may be reduced in three-body decays, and experimental tests of this 
hypothesis - is examined in Sec. 5. We conclude in Sec. 6 . 

2 . Dalitz Plots 

Three-body B decays are usnally described nsing a Dalitz plot. Consider the decay 
B —)■ P 1 P 2 P 3 , in which each pseudoscalar Pi has momenta pi. One can constrnct 
the three Mandelstam variables Sij = {pi +pj)^, where pi is the momentnm of each 
Pi. These are not independent, bnt obey <§12 + -S 13 4- S 23 = + ml + ni^ + m|. 

The B —)■ P 1 P 2 P 3 Dalitz plot is a measnre of the decay rate as a fnnction of two 
Mandelstam variables. 

In the present paper we focus on the decays B^g —)■ Rs{pi)h~^{P 2 )h~{p 3 ) {h = 
R,7t). At the qnark level, the hnal states / = Rst^^t^~ and RsR~^R~ are self- 
conjngate. However, when the momenta are considered, one has f ^ f. The point 
is that the CP conjugate of / = Rs{Pi)h^{P 2 )h~{p 3 ) is / = Rs{Pi)h~{p 2 )h+{p 3 ), 
where pi is pi with the direction of the three-momentnm reversed. Note that reversing 
the direction of the three momenta does not affect the Mandelstam variables, since 
^ij = {Pi PPj)^ = {pi +Pi)^ = ^ij- Thus, in this case the difference between / and / 
arises from an exchange of the indices 2 and 3. 
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The distinction between / and / must be kept in mind throughout the paper. 
Because / is self-conjugate at the quark level, both and can decay to it, and 
similarly for /. Now, at different points in the analysis we consider the direct CP 
asymmetry. However, because / 7 ^ /, there are two of these. One compares —)■ / 
and B^ ^ f decays, the other B^ ^ f and B^ ^ f. Things are similar for the 
indirect CP asymmetry, which arises because both B^ and B^ can decay to the same 
hnal state. Thus, one indirect asymmetry involves the interference of the amplitudes 
for B^ ^ f and B^ —)■ /, while the other involves the interference of A{B^ f) and 
A{B^ ^ /). 

3. U-Spin Relation 

In this section we discuss the U-spin relation that is central to our method for ex¬ 
tracting 7 . We begin by reviewing the relation for two-body decays. 

3.1 Tw^o-body decays 

Consider a pair of H —)■ PP decays whose amplitudes are related by U-spin reflection 
(d -H- s). (This discussion follows Ref. [28].) One is a 6 —)■ d decay, the other 6 —)■ s. 
There are hve such pairs [16]: (R° —)■ n^n~, —)■ K^K~), (R° —)■ 7r~^K~, B^ —)■ 

7r-i^+), (R+ ^ K+K^, B+ 7r+iP°), (5° ^ ^ (5° ^ K+R-, 

B^ —)■ 7r’''7r“). 

The b ^ d amplitude can be written 

A, = AuV:,Vud + A,v:,v,d + AtV:,Vu 

= - A,)v:,v^d + (die - A)UetKd 

= v:^VudTd+v:^v,dPd. ( 3 . 1 ) 

In the above, the Ai each represent a linear combination of diagrams, and we have 
used the unitarity of the CKM matrix (y*^Vud + Vcb^cd + Vth^td = 0) to write the 
second line. and Pd are simply the quantities that are multiplied by the given CKM 
matrix elements - they do not represent individual “tree” and “penguin” diagrams. 
The 6 —)■ s amplitude can be written similarly: 

A = + V:,V,sPs . (3.2) 

The CP-conjugate amplitudes Ad and As are obtained from the above by chang¬ 
ing the signs of the weak phases: 

Ad = + K,V;,Pd , A = KbK^Ts + K,V*Ps . (3.3) 

We then have 

l^r - lAr = 4Im(lC,14,K*l/aim(T,P,*) , 

I Ap - 77 = 4 Im(V:,K.K,V;j Im(T.P,-) . (3.4) 
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Now, the unitarity of the CKM matrix implies [29] 

, (3.5) 

and in the U-spin limit we have 


T, = T. 


Ph = P. 


(3.6) 


U-spin symmetry therefore leads to a relation between the b ^ d and 6 —?• s decays: 


\Ad\^-\A,\^ = -[\As\^-\As\^] . (3.7) 


In general, there are four observables in the b ^ d and b ^ s processes: the 
branching ratios Bd and Bg, and the direct CP asymmetries and . Eq. ( |3.7]) 
implies that these are not independent, but obey [9,28] 




(3.8) 


Thus, there are only three independent observables. 


3.2 ^ —>■ Ksh^h decays 

We now turn to B^^ —)■ Ksh^h~ decays. For definitiveness, we focus on the pair 
(Bg —)■ iP^TT+TT” (b ^ d), B^ ^ KsK^K~ (b —)■ s)), but the results can be equally 
applied to (5° —)• KsK^K~ (b ^ d), B^ ^ (b —)■ s)). 

As discussed in Sec. ||, one must pay attention to the momenta of the final-state 
particles. Let us define fd = Ks(pi)n+(p2)7i~(p3) and fd = Ks(Pi)'^^(P 3 )'^~(P 2 ), and 
similarly for fg and fg. Now consider Ad = A(Bg —)■ fd) and Ag = A(B^ —)■ fg). The 
decay amplitudes Ad and Ag are again given by Eqs. (|3.1|) and (|3.2|), respectively, 
and are repeated below for convenience: 


= V:,VudTd + VXV^dPc 


cb^cd^d , 


A = V:,V^gTg + VXV^gPg 


cb csA S • 


(3.9) 


As these are three-body decays, Td^g and Pd^s are all momentum-dependent. This 
means that Td takes different values at different points of the Dalitz plot, and similarly 
for Tg and Pd,s. For the CP-conjugate amplitudes, we have 

.Arf = VubVfdfd + VcbVfdPd , A = VubVfJg + V,bV,lPg . (3.10) 


Because the final states in the CP-conjugate decays are not the same as in the decays 
(p 2 and p 3 are exchanged), Td ^ Td, and similarly for Tg and Pd^s- 
We then have [30] 

= 2 lm(V:,VudVcbV:i) lm(TdP: + ffPd) , 
lAP - lAP = 2lm{V:,VusV,,V:,)lm{TgP: TTfPg) . (3.11) 
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In the U-spin limit we have Td = Tg, Pd = Ps, = Ts, Pd = Ps, and the U-spin 
relation of Eq. (|3.7|) is reproduced. However, since the amplitudes themselves are 
now momentum dependent, this relation holds at each point in the Dalitz plots. 

As in the two-body case, the U-spin relation implies a relation among the ob¬ 
servables, similar to Eq. (|3.8| ). This relation involves —)■ / and —)■ / decays, 

and can be written as 


T{K)h 


,CP 


Here, and bg are, respectively, the direct CP asymmetry and branching ratio 
dehned locally, i.e., at a particular Dalitz-plot point. They are both momentum- 
dependent quantities. 

The analysis can be repeated for the case where Ad = A^B'^ —>■ fd) and As = 
A{B^ —>■ fs). Here we have 


Ad = v:,v^dfd+v:,v,dPd , 

and 

Ad = Vubv:dTd+Vcbv:dPd , 


A = v:,v^sTs + v:,v,sPs . 

(3.13) 

A = KfeCT, + v,,v:sPs . 

(3.14) 


Once again, the U-spin relation of Eq. (^.7| ) is reproduced. And there is a relation 
like Eq. ( p.l2| ) among the observables. This relation involves 5° —)• / and 5° —)• / 
decays. 

The point here is that, for three-body decays, there are two U-spin relations 
among the observables. These involve the same momentum-dependent hadronic pa¬ 
rameters. 


4. Extraction of 7 

Here we present the details of how 7 can be extracted from a U-spin analysis of 
B^^^ —)■ Ksh'^h~ decays. We begin with a review of the method for two-body decays. 

4.1 Two-body decays 

The method proposed by Fleischer for extracting 7 from 5° —)■ K^K~ and B^ —>■ 
TT+TT” [9] works as follows. The amplitude for the b ^ d decay {B^ —)■ tt+tt”) is given 
in Eq. (p.l| ), which can be written 

- \V:,V,d\Pd , (4.1) 

where we have used \Vcd\ = —Wd- The amplitude for the b ^ s decay (i?° —)■ K^K~) 
can be written similarly: 




(4.2) 
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In the U-spin limit, we have = Tg = T and Pd = Pg = P■ Assuming that the 
magnitudes of the CKM matrix elements are known, Ad and Ag each contain the 
same four unknown parameters: |T|, \P\, their relative strong phase, and 7 . 

Above [Eq. ( p. 8 |) ], it was noted that the branching ratios and the direct CP 
asymmetries of these two decays are not independent. Thus, 7 cannot be extracted 
from the measurements of these observables alone, since there are more unknown 
theoretical parameters (four) than observables (three). However, if the indirect CP 
asymmetries in both —)■ 7 r’'' 7 r“ and —)■ K~^K~ are also measured, and values 


for the B^-B^ and B^-B^ mixing phases (/? and /3s, respectively) are taken from 
independent measurements, there will be more observables (hve) than unknowns, 
which will allow 7 to be extracted. 


4.2 ^ —>■ Ksh^h decays 

A similar logic can be applied to three-body decays. However, care must be taken 
in identifying the observables to be used, and in establishing how these observables 
depend on the unknown theoretical parameters. 

The point is the following. If a hnal state / is self-conjugate at the quark level, 
both B^ and B^ can decay to it. In the case of two-body decays, the fact that / 
is self-conjugate implies that / = /, so that the two decays 5°, 5° —)■ / must be 
considered. However, as noted in Sec. ||, for three-body decays, a self-conjugate / 
still has / 7 ^ /, since / and / correspond to different points of the Dalitz plot. In this 
case, the analysis must consider the four decays 5°, B^ —)■ /, /. The time dependence 
of two-body decays has been analyzed in Refs. [31,32]. Below we adapt this analysis 
to three-body decays. 

In the presence of B^-B^ mixing, the B^ and Bh states (L is light, H is heavy) 
are mixtures of B^ and B . The physical time-dependent neutral R-meson states 
can then be expressed as 




R 


B 


phys 


it)) = 7-W 7")+ /+(<) I®") 


(4,3) 


Here 5phys(^) (^hys(^)) is the state that is a B^ (i?°) at f = 0. In the above, 
q/p = where (j)M is the weak phase of the mixing (the B^-B^ and 

mixing phases are (3 and /3s, respectively), and 


f+{t) = cos (A/if/2) , /_(f) = e-'(”^-'^/'i*isin(A/if/2) , (4.4) 


with 


m = {niH + mL) /2 , Am = mn — mi , 
r = {rH + rL)/2 , Ar = rH-rL, 
A/i = Am — iAT /2 . 


(4.5) 






The decay amplitudes are then given by 


where 


(/|Bph,.(i)> = (/|B"> (/+(*)+A/-(*)) . 

(7|b*.(*)> = ^ (7|5”) (f+m + 

(/iKi^ys(t)} = ^ (/|B“> (/-(*) + A/+(«)) . 

(jlKiyM) = (71^) (/-(*)A+ /+(*)) . 



'/|b“ 


g p 

A = -X , A = -X 

P Q 


(4.6) 


(4.7) 


In Ref. [31] the assumption is made that AT = 0. In Ref. [32] it is noted that 
AT is nonzero in decays. Our expressions below therefore allow for a nonzero 

AT. 

The decay rates are proportional to the squares of the amplitudes, which take 
the form 


^ f) 

|-MP(5physW ^ /) 
\M\\Bl^^^{t) ^1) 
|-MP(5physW^7) 


|Afe (1 — \x\^) cos(Amt) + (1 + \x\^) cosh(Art/2) 


-2 Im(A) sin(Amt) + 2 Re(A) sinh(Art/2) 




2-Vt 


-(1 — |a;|^) cos(Amt) + (1 + |a;p) cosh(Art/2) 


+2 Im(A) sin(Amt) + 2 Re(A) sinh(Art/2) 




|A|^e —(1 — |a;|"') cos(Amt) + (1 + |a;|"') cosh(Art/2) 


+2 Im(A) sin(Amt) + 2 Re(A) sinh(Art/2) 


-\Afe ^ (1 — |a;f) cos(Amt) + (1 + |a;f) cosh(Art/2) 


-2 Im(A) sin(Amt) + 2 Re(A) sinh(Art/2) , (4.8) 


where A = (/|R°), A = (^f\B^^, and we have used \q/p\ = 1. 

With the squares of the amplitudes in hand, we can now obtain expressions 
for the observables. Before doing so, there is one point that must be mentioned. 
Although we have referred to measurements at different points of the Dalitz plot, 
in practice it is only possible to make measurements in bins, i.e., over areas of the 
Dalitz plot centred at different points. The observables will then involve integrals 
over the Mandelstam variables representing these bins. 



Using the first two equations of Eq. (|4.8|) , we can now construct the time- 
dependent CP-averaged rate and the CP asymmetry for the final state /: 


m = ^ /) + ^ /)), 

dsi 2 ds 23 [(l -I- |a;p) cosh(APt/2) 


bin 


2Re(A) sinh(APt/2))] , 


(4.9) 


^cpif) — 


r(Bghy,(i) ^ /) - nBluyM ^ /) 

r(B”h„(*) ^ /) + n'Kl.y.W ^ /) ’ 

Jf dsudS 23 |4l|^ [(1 — |a:p) cos(Ami) — 2Im(A) sin(Ami)] 

bin 

// dsi 2 ds 23 \ A\^ [(1 -|- |xP) cosh(APt/2) -|- 2Re(A) sinh(APt/2)] 

bin 


.(4.10) 


In P(f), one does not distinguish 5phys(^) decays, whereas one does in 

Ac pit). Thus, as usual, the measurement of the CP asymmetry requires tagging. 


A comment should be made about Eq. ( 4.10 ). The direct CP asymmetry com¬ 
pares 5° —;• / and —)■ / decays. Because / = / in two-body decays, there 
one refers to the coefficient of cos (Amt) as the direct CP asymmetry. However, in 
three-body decays, because f ^ f, the situation is different. Here the coefficient of 
cos (Amt) compares ^ f and B^ ^ f decays, and so it is not actually a CP 
asymmetry. 

In the above definitions there appear to be four observables, namely the coeffi¬ 
cients of cos(Amf), cosh(APf/2), sin(Amf), and sinh(APf/ 2 ), as can be determined 
from P(f) and the numerator of Acp{t). However, these coefficients are not all inde¬ 
pendent, as can be seen in the following identity: 


\AWi + \x\^) -\AWi - \x\^) = 2\A\^\x\^ = 2|AnAp 


= 2|A|2(Re(A)2 + Im(A)2) 


(4.11) 


There are therefore only three independent observables. 

One can perform a similar analysis using the last two equations of Eq. ([4.8|) . In 
this way one constructs the time-dependent CP-averaged rate and the CP asymmetry 
for the final state /. There are again three independent observables. Thus, for a 
given —)■ Ksh^h~ decay, there are a total of six observables: three each for the 

hnal states / and /. 

We now turn to the question of the number of unknown theoretical parameters, 
focusing on the decay pair —)■ iP^Tr+vr” {h —?• d) and B^ —>■ KsK^K~ (6 —)■ s). 

Consider Erst the b ^ d decay. The amplitudes for the various B^,B^ —)■ /, / 
decays are given in Eqs. (|3.9|) , (|3.10|) , (|3.13|) and ( p.l4|) . There are eight unknown 
parameters: \Td\, \Pd\, \Td\, \Pd\, their three relative strong phases, and 7 . With only 
six observables, 7 cannot be extracted. 


9 













Td = T, = T and = P. 

eight unknown parameters. (As before, it is assumed that the 


This can be remedied by also considering the U-spin conjugate 6 —)■ s decay 
—)■ KsK^K~. Its 5°, 5° —t /, / amplitudes are also given in Eqs. (|3.10|) , 

( p.l3|) and (p.l4|) . Here too there are eight unknown parameters: \Ts\, |Ps|, \Ts\, 
|Fs|, their three relative strong phases, and 7 . However, in the U-spin limit, we have 

Thus, the two decays are described by the same 

° B^ and 

mixing phases are taken from independent measurements.) But there are now twelve 
observables, six for each of B^ —)■ Ks7r~^'n'~ and B^ —)■ KsK^K~. On the other hand, 
it was noted in Sec. I 


I that there are two U-spin relation among the branching ratios 
and direct CP asymmetries of the 6 —)■ d and 6 —)■ s decays. Still, this leaves ten 
independent observables, which is more than the number of unknown parameters. 
Thus, assuming again that the magnitudes of the CKM matrix elements are known, 
7 can be extracted. 

It must be mentioned that this method introduces a new systematic error. We 
have argued above that since the number of observables is greater than the number of 
unknowns, 7 can be extracted. But this only works if all the observables are functions 
of the same unknowns. And because the measurements must be made using bins of 

n=:0' 


the Dalitz plot, this does not hold exactly. Writing A x = {f\B ) = A, from Eqs. 
and ( |4.10|) we have 


BR oc JJ dsi 2 ds 23 (|Ap -|- |Ap) , 

bin 

^dir ^ jj dSi2dS23 (|A|^ - lAp) , 


oc Jj dsi2ds23lm[{q/p)A*A] . (4.12) 


If we dehne 

bin 



JJ dsi2ds23 Ap 

bin 

= |Af , 

J J dsi2ds23 Ap = A'p , 

bin 

(4.13) 


we see that both BR and are functions of A' and A'. However, is not. 

We must make the approximation that 


jj dsi2ds23lra[{q/p)A* A] ~ Im[(g/p)A'*A'] , (4.14) 

bin 

and this introduces a systematic error. The above holds exactly for a single point 
of the Dalitz plot. Thus, the smaller the bins are, the better is the approximation, 
leading to a smaller systematic error. On the other hand, smaller bins lead to larger 
statistical errors. The bin size must therefore be chosen to minimize the total error. 
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5. U-Spin Breaking 

As noted earlier, the method of combining measurements of decays related by U 
spin to extract 7 was originally proposed in the context of two-body decays [9]. 
Here, there is a theoretical error due to unknown U-spin-breaking effects. This same 
difficulty arises when applying the method to three-body decays. In this section we 
examine the question of U-spin breaking as pertains to three-body decays. 

The method described in the previous section for extracting 7 involves combining 
measurements of pairs of three-body decays related by U-spin, such as —)■ Ks7r~^7r~ 
(b —)■ d) and —)■ KsK^K~ (6 —)■ s). This method applies at a particular pair 

of Dalitz-plot points (bins). By repeating this analysis for all points, this provides 
multiple measurements of 7 . These can then be averaged over the entire Dalitz plot, 
reducing the statistical error. 

In the presence of U-spin breaking, the extracted value of 7 , 'jext, will differ from 
its true value, 'jtrue- Now, there are several different U-spin-breaking parameters. 
However, these parameters are all momentum dependent. Thus, their effect on the 
extracted value of 7 will vary from point to point on the Dalitz plot. That is, if 

dexf 'Itrue -N , (5.1) 


it is likely that iV > 0 at some points, and < 0 at others. In this case, averaging 
over all Dalitz-plot points will also reduce the effect of U-spin breaking, so that 
i'yext)avg will approach 'ytrue- If this occurs, the main theoretical error of the method 
will be signihcantly reduced. 

Still, while this is a nice idea, how can we be certain that it is happening? 
Fortunately, there is a way of experimentally testing whether or not this behaviour 
is present in three-body decays. In Eq. ( 3.12|) it was shown that there is a relation 
among the observables of two decays related by U-spin reflection. Writing 


.CP 


.CP 


r{Bj)bs 

rmbd 


1 = n' , 


(5.2) 


we have n' = 0 in the U-spin limit. By measuring 6^,5 and , and constructing the 
above ratio at each Dalitz-plot point, it is possible to experimentally determine if an 
average over all points leads to n' —)■ 0 . 

The above test requires a Dalitz analysis. A simpler test of U-spin breaking can 
be obtained by separately integrating the numerator and denominator of Eq. ( 0 ) 
over the kinematically-allowed regions of the Dalitz plots (denoted by DP): 


II dsi2ds23a^^b 

DP_ 

II dsi2ds23a^^b 


S 


d 


- 1 


AJ'’ T{B«)Ba 


(5.3) 
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Unlike n', which is dehned using momentum-dependent quantities, N' depends only 
on integrated quantities, and hence does not depend on hnal-state momenta. Once 
again, we have = 0 in the U-spin limit. 

The above tests can be carried out using the measurements of —)■ Ksh'^h~ 

decays. However, it is not necessary to wait until these are made. Other pairs of 
three-body decays related by U spin are (i) 'k^K^K~ (b —)■ d) and 

(b ^ s), and (ii) tt+tt+tt" (b —)■ d) and B^ —)■ K^K^K~ (b ^ s). 

In Ref. [33], group theory is used to write the factor n' of Eq. (|5.2|) for these decay 
pairs in terms of U-spin-breaking parameters. These parameters take into account 
all U-spin-breaking effects, such as differences in the masses of the tt and K mesons, 
differences in the properties of the resonances contributing to the decays (e.g., p, 
0), etc. It is found that, to hrst order, n' is proportional to a linear combination 
of such parameters. That is, in the presence of U-spin breaking, n' 7 ^ 0. However, 
the U-spin-breaking parameters are momentum-dependent. Using the same logic as 
before, it would not be surprising to hnd n' > 0 at some points and n' < 0 at others. 
If so, the average over all Dalitz-plot points will reduce the effect of U-spin breaking 
in the above relation. 

These B^ decays have recently been measured by LHCb [34,35]. In Ref. [30], 
the U-spin relation of Eq. ( |5.3|) is tested using data integrated over the Dalitz plot. 
We have updated these results with more recent data from Refs. [36]. The updated 
results are shown in Table |^. Unfortunately, at present the results are simply not 
precise enough to draw any conclusions. When the data improve, we will have a 
better idea of whether averaging (or integrating) over the Dalitz plot reduces the 
effect of U-spin breaking. 


Asymmetry ratio 

U-spin 

prediction 

LHCb 

result 

^ 'k+K+K-)/A^P{B+ AT+vr+TT-) 
A^p\b+ 7r+7r+7r-)/A^^(R+ ^ K+K+R-) 

-10.2 ± 1.5 

- 2.2 ± 0.2 

-4.9 ±2.0 

-1.6 ±0.5 


Table 1: U-spin predictions for asymmetry ratios compared with LHCb measurements. 


Finally, another source of U-spin breaking arises from the fact that vr^ and 
do not have the same mass, and similarly for B^ and i?°. This results in a difference 
between the kinematically-allowed phase space for a decay and that for its U-spin 
partner. Due to this difference, there will be regions of the Dalitz plots where the 
observables dehned in Sec. ^ can be obtained only for one of the two decays being 
compared. These regions must be excluded from the analysis, since our method for 
extracting 7 works only for those regions of the Dalitz plots where the two decays 
have overlapping kinematically-allowed regions. 
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6. Conclusions 


In 1999, R. Fleischer proposed a method for extracting 7 using a pair of two-body 
decays whose amplitudes are related by U-spin symmetry (d <H- s) [9]. It involves 
combining the measurements of the branching ratios and CP asymmetries, both 
direct and indirect, of the two decays. These decay amplitudes include penguin 
diagrams, which may receive important (loop-level) contributions from new physics. 
If so, the value of 7 extracted using this method will disagree with its current value, 
which is obtained using tree-level decays. 

In the present paper we adapt this method to charmless B —)■ PPP decays (P 
is a pseudoscalar meson), specihcally —)■ Ksh~^h~ {h = K,n). Time-dependent 

Dalitz analyses of the three-body decays can be used to measure the branching 
fractions and CP asymmetries. Note that it is not necessary to perform an isobar 
amplitude analysis of the Dalitz plot. We show that there are more observables 
than unknown theoretical parameters, so that 7 can be extracted by htting to the 
observables. The decay amplitudes for three-body decays depend on the momenta 
of the hnal-state particles. The method applies to each point of the Dalitz plot, and 
thus constitutes many independent measurements of 7 . 

The main source of theoretical error in the extraction of 7 , which also applies 
to the method with two-body decays, is U-spin breaking. However, three-body de¬ 
cays offer the potential to reduce this error. The U-spin-breaking effects are also 
momentum-dependent. As such, the difference between the extracted value of 7 and 
its true value may well vary, in both magnitude and sign, from point to point in the 
Dalitz plot. If this is the case, then averaging over the Dalitz plot will reduce the 
error due to U-spin breaking. 

It is possible to test experimentally whether or not this behaviour is present 
in three-body decays. In the U-spin limit, there is a relation among the branching 
ratios and direct CP asymmetries of the two decays that are related by U spin. This 
applies to the decays —)■ (b d) and B^ —)■ A'+7r’''7r“ (b ^ s), and 

B'^ —)■ 7 r''' 7 r''' 7 r“ (b d) and B^ K^K^K~ (b s), all of which have been 
measured. Unfortunately, the current data on these decays still has large errors, so 
that it is unclear whether U-spin breaking is small when averaged over the Dalitz 
plot. Future precision data in these channels will be able to clearly show the size of 
this U-spin breaking. 
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